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STRESS CONCENTRATIONS IN FILAMENT-STIFFENED 
SHEETS OF FINITE LENGTH 

By W. B. Fichter 
Langley Research Center 

SUMMARY 

A simple model of filamentary composite material is employed to investigate stress 
concentrations in a filament- stiffened sheet of finite length. The model is composed of a 
single layer of parallel, tension- carrying filaments embedded in a shear- carrying matrix. 
The sheet is of finite length in the filament direction and of infinite length normal to the 
filament direction. Filament stress- concentration factors are calculated as functions of 
the number of broken filaments and a length- stiffness parameter for the cases of uniform 
normal edge load and uniform normal edge displacement. In the uniform-edge-load case, 
the stress- concentration factors are found to increase with decreasing filament length. 

The opposite effect is noted in the uniform-edge-displacement case where, in addition, the 
stress-concentration factor is found to have an upper limit which is fixed by the value of 
the length- stiffness parameter. 


INTRODUCTION 

In the design and application of aerospace structures the problem of stress concen- 
tration is one of continuing importance. In recent years conventional engineering metals 
have been replaced in some weight- critical applications by filamentary composite mate- 
rials. These materials behave somewhat differently from conventional metals when 
weakened by flaws and, hence, different analytical models must be used to investigate 
their behavior. 

One simplified model for filamentary-composite-material behavior which has been 
found useful is the "stringer-sheet" model employed in reference 1 for investigating sta- 
tic and dynamic stress concentrations around broken filaments. This model consists of a 
sheet of parallel, equally spaced, tension- carrying filaments embedded in a shear- 
carrying matrix. In references 2 and 3 the model has been used to investigate the inter- 
action of neighboring groups of broken filaments, the matrix shear forces, and the 
recovery of load by broken filaments. In reference 4 a similar model has been employed 
to study stress concentration in two-dimensional arrays of equally spaced filaments. 



In all of these studies attention has been restricted to bodies of infinite extent. It 
is well known, however, that results for infinite bodies cannot always be applied indis- 
criminately to similar problems involving finite specimens. It is desirable, therefore, to 
assess the influence of specimen size on stress concentration. In the present paper the 
stringer- sheet model is employed to investigate static stress concentration due to broken 
filaments in a filament- stiffened sheet of infinite length normal to the filament direction 
and of finite length in the filament direction. Stress-concentration factors are found for 
a central cut across various numbers of adjacent filaments and for a wide range of values 
of a nondimens ional length- stiffness parameter. 

SYMBOLS 


d 

EA 

Gh 

K r 

k 


l 

M 


M. 


n 


m,n,r 


N n 

p n 

P 


Pn 


r 


filament spacing 

extensional stiffness of a filament 
matrix shear stiffness per unit length 
stress- concentration factor for r broken filaments 
length- stiffness parameter, l 
half-length of filament sheet 
transformed influence functions 

displacement of nth filament for influence- function solution 
indexes 

load in nth filament for influence- function solution 
dimensionless load in nth filament, p n /p 
filament load parameter 
load in nth filament 
number of broken filaments 
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FA 

U n dimensionless displacement of nth filament, 

pi 

u n displacement of nth filament 

x coordinate parallel to filament 

rj dimensionless coordinate parallel to filaments, 

6 transform variable 

Superscript: 

00 denotes values for infinite- length sheet 


u 


n 


x 

l 


ANALYSIS 


The finite-length sheet under investigation is shown in figure 1 along with the coor- 
dinate system and some notation. The filaments are separated by a constant distance d 
and are numbered from - °° to °° from left to right. The coordinate along the fila- 
ments is x and the displacement of the nth filament at location x is denoted by u n (x). 
The force in the nth filament, taken to be positive in tension, is p n (x) and is given in 
terms of u n by 


p n = EA 


dun 

dx 


( 1 ) 


where EA is the extensional stiffness of the filament. The matrix shear force per unit 

r i_ 

length between the nth and (n + l)th filaments is given by — (u n+ j - u n j, where Gh 
is the shear stiffness of the matrix. Static equilibrium of an element of the nth filament 
requires that 

2 

E A TT + IT ( Un+1 " 2un + u n- 1) = 0 (2) 

QX 

In figure 1 filaments 0 and 1 are shown broken at x = 0. In general, for r broken 
filaments, let 0 S n s r - 1 denote the broken filaments. Then, the appropriate bound- 
ary conditions along the line containing the cut (x = 0) are 
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p n (0) = 0 


(0 £ n £ r - 1) 


u n (0) = 0 


(n < 0 or n £ r) 


( 3 ) 


Two sets of boundary conditions at the edges of the filamentary sheet are of interest. 
These are as follows: (I) uniform normal displacement, and (II) uniform normal load. 
In terms of the pertinent variables, the two conditions are 


Problem I: u n (±Z) 



(4) 


and 


Problem II: p n (±Z) = P (5) 

(Note that these two conditions would be equivalent if no filaments were broken.) Prob- 
lems I and II are defined by equations (2) and (3) along with equations (4) and (5), 
respectively. 


Nondimensionalization 

In order to cast the problems in a more convenient form, let 


Pn = PP n 




u n - u n > 


x =lri 


( 6 ) 




Then, equation (2) becomes 


d 2 U r 


dry* 


?+k 2 (U 1M . 1 -2U n + U n _ 1 ) =0 


(7) 


where k is a length- stiffness parameter defined by l For the two problems 

considered here the complete sets of boundary conditions become 

P n (0) = 0 (Ognsr-l) 

U n (0) =0 (n < 0 or n s r) 


( 8 ) 
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and 


Problem I: U n (±l) = ±1 (9a) 

and 

Problem II: P n (±l) = 1 (9b) 

The dimensionless forces and displacements are related by 
dU n 

Pn(?7) = “ST (r?) (10) 

Solution of problems I and II is complicated by the fact that the boundary conditions 
at x = 0 are mixed; that is, load is prescribed over one part of the boundary, with dis- 
placement prescribed over the remainder. To overcome this difficulty the influence- 
function technique introduced in reference 1, and later used in references 2 and 3, is 
again employed. It is assumed that the force and displacement in the nth filament have the 
form 


and 


GO 

P n (r?) = 1 + £ N n _ m (77)U m (0) 


oo 

U n (T7)=7]+ ^ M n _ m (i7)U m (0) 


m= 


where 


dM n (77) 

dr] 


= N n (77) 


The quantities M n and N n are, respectively, the displacement and load in the nth fila- 
ment due to a unit displacement of the zeroeth filament at rj = 0. Application of the bound- 
ary conditions at 77 = 0 (eq. (8)) yields 
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r -1 

P n (rj) = 1 + £ N n _ m (rj)U m (0) 
m =0 

r -1 

U n (77) =V+ ^ M n _ m (7 1 )U m (0) 
m=0 

since U m ( 0 ) = 0 for other values of m and 
r -1 

° = 1 + £ N n _ m (0)U m (0) (Osngr-1) (12) 

m=0 

which expresses the condition of zero load at the breaks in the r filaments identified by 
n = 0 to n = r - 1. It should be noted that because of symmetry only the upper half of 
the sheet (0 s -q = 1) need be considered. Equations (12) constitute a set of r linear 
alegebraic equations for the r unknown U m (0). They can be solved and their solution 
substituted into equations (11) to yield the entire solution. Before this can be done, how- 
ever, the influence functions for the particular problem of interest must be determined. 
Because of the similarities between the influence functions for the two problems, they are 
determined here simultaneously in order to avoid unnecessary repetition. 

Determination of the Influence Functions 

The influence functions M n and N n are, respectively, the nondimensional dis- 
placement and load in the nth filament. These result from displacing the zeroeth filament 
a unit amount at rj = 0 while maintaining zero displacement at r) = 0 for all other fila- 
ments, and also while maintaining at the edge rj = 1 either zero displacement (prob- 
lem I) or zero load (problem II) in all filaments. 

For 17 £ 0 the problem can be stated as 
2 

^sr + k2 ( M n +1 " 2M n + M n-l) = 0 < 13) 

dr^ v ' 

with the conditions 
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M n (0) = 1 
M n (0) = 0 


(n = 0) 
(n* 0) 




M n (l) 

dM 


( 14 ) 


= 0 


dry 


2(1) 




Here and in the remainder of this report the brace symbol ^ } means that the upper term 
corresponds to problem I (uniform normal edge displacement) and the lower term corre- 
sponds to problem II (uniform normal edge load). 

Let 


M (ry,e) = ^ M n ( 7 y)e me 

n= 

or, inversely, 

M n (i7) = i- C m-n,6)e ind de 

2tt J-v 

Application of equation (15) to equation (13) yields 

- 4k2 sin^ ~ M = 0 
^2 2 

Similar treatment of the boundary conditions in equations (14) yields 
A 

M(O,0) = 1 


M(M) 


If 1 »•*> 


> = ° 


The solutions satisfying equations (17) and (18) are found in 


M(ry,0) = cosh (2k?y sin ^ 


- 


coth (2k I sin £ 


. 6 
sm — 
2 


tanh ^2k 


sinh ( 2k?; 


Sin t 


(15) 


(16) 


(17) 


(18) 
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Substitition of these solutions into equation (16) gives 




and differentiation yields 
n r 


p 71 




r / 


\ 

~ 


/ 

sin 

2 

\ , 

coth ( 2 k 

sinf' 

/ 

sin 1 1) 

i 

cosh [ 2 k T) 

)- 

V 

tanh [ 2 k 

2 y 

sin l) 

Ssinh ( 2 k 77 


L 





L 

J 



e ln9 de 


N„(u) = | 

71 


-7T 


sinh [2k rj 


. 0 
sm -r 
2 




coth ^2k 
tanh /2k 


. 0 
sm 2 

sin 

2 


\cosh 2k?] 


sin | 




sin | 


e inS de 


Finally, making an appropriate change of variables and recognizing the symmetry proper- 
ties of the integrand yield 


au pir/2 [coth (2k sin 0) 

N n (0) = - — \ cos 2n0 sin 0 < > d0 

v J 0 (tanh (2k sin 0)1 


(19) 


Stress-Concentration Factors 

The maximum force occurs at 77 = 0 in the first intact filament adjacent to the 
broken ones. Therefore, for a cut across r filaments the stress- concentration factor 
K r is given by P r (0) which is equal to P_j(0) because of symmetry. Thus, 

r-1 

K r s P r (0) = 1 + £ N r _ m (0)U m (0) (20) 

m=0 

where values of N n (0) for problems I and II are obtained from equations (19). The 
values of U m (0) are obtained from equations (12) after substitution of the appropriate 
values for N n (0). 

Closed-form evaluation of the influence functions N n (0) given by equations (19) 
does not appear feasible, although asymptotic expressions valid for small values of n 
and either very small or very large values of k are not too difficult to obtain. There- 
fore, the influence functions have been evaluated by numerical integration on a digital 
computer, with the asymptotic expressions being employed in a few cases merely as a 
check on the numerical results. In addition, the solution of equations (12) and calcula- 
tion of the stress- concentration factors from equation (20) have been included in the auto- 
matic routine. 
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RESULTS AND DISCUSSION 


For both of the problems considered, stress-concentration factors have been calcu- 
lated for various numbers of broken filaments ranging from 1 to 100 and for a wide range 
of values of k, the length- stiffness parameter. The results for problem I are presented 
and discussed first. 


Problem I - Uniform Normal Edge Displacement 

The stress-concentration factor K r due to various numbers of broken filaments (r) 
is plotted in figure 2 as a function of k, the length- stiffness parameter. The results 
for small numbers of broken filaments are shown in figure 2(a), and for larger numbers in 
figure 2(b). Also shown are the stress- concentration factors for the infinite -length sheet 
(see ref. 1) to which the present finite-length results are asymptotic from below for 
large values of k. 

For very small values of k, it is seen in figure 2(a) that the stress- concentration 
factor is only slightly greater than one and is independent of the number of broken fila- 
ments. As k increases, however, the number of broken filaments becomes increasingly 
important. For example, for k = 0. 5, no effect of the number of broken filaments is 
evident; whereas for k = 10, different numbers of broken filaments yield significantly 
different stress -concentration factors, as long as the number of broken filaments is less 
than about 20. For even larger values of k, the number of broken filaments is of signif- 
icance over a still wider range. This range is still limited, however, for any finite 
value of k. 

In an experimental program of tensile tests of filamentary composites with cuts, it 
is likely to be k, the length- stiffness parameter, rather than r, the number of cut fila- 
ments, which would be held constant throughout a number of tests. Therefore, it would be 
informative to present the results in a different way. This is done in figure 3, where the 
stress- concentration factor is plotted as a function of the number of broken filaments, 
for various values of k. Although the results are valid only for discrete values of r, 
they are plotted as continuous curves for ease of illustration. It can be seen that for a 
given value of k the stress-concentration factor ultimately attains a value which cannot 
be exceeded, no matter how many filaments are broken. For example, for k = 4, the 
stress- concentration factor has an upper limit of about 2.10; and for k = 10, an upper 
limit of about 3.25. Wherever feasible it would appear to be desirable in many applica- 
tions to design a composite structural component in such a way that an acceptable upper 
limit could be imposed on the filament stress- concentration factor. This is analogous 
to the "fail-safe" philosophy of aircraft design. Within the limitations of the model, this 
upper limit could be imposed by prescribing a value of k and then suitably adjusting 
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the relative values of length, matrix shear stiffness, filament extensional stiffness, and 
filament spacing. 


Problem II - Uniform Normal Edge Load 

Figure 4(a) contains results for small numbers of broken filaments; the results 
for larger numbers are presented in figure 4(b). In contrast to the results for problem I, 
the stress- concentration factors for the uniform- edge -load problem are greater than 
those for the infinitely long sheet, and they approach them from above as k becomes 
large. For very small values of k, it can be shown that the stress -concentration factor 
K r approaches the value 1 + This rather severe stress concentration is largely of 
academic interest, however, since uniformly loaded edges are rarely encountered in 
practice. In the event of their occurrence, however, the present results suggest that 
modifications should be made, perhaps by reinforcing the edges of the sheet, in order 
that the boundary conditions would correspond more closely to those of problem I. 

CONCLUDING REMARKS 

Filament stress-concentration factors have been calculated for two problems con- 
cerned with the stretching of a finite-length sheet of filamentary material weakened by 
the presence of a group of broken filaments. For the problem of uniformly loaded edges, 
the effect of finite length is manifested by an increase in the stress-concentration factor 
over that for the infinite -length sheet. However, in the problem of uniform edge dis- 
placement, the stress-concentration factor is lower for the finite-length sheet than for 
the infinite -length sheet. In addition, it is seen that to each given value of the length- 
stiffness parameter there corresponds an upper limit which the stress- concentration 
factor cannot exceed, regardless of the number of broken filaments. This limit should 
be of interest to designers working within the constraints of a maximum allowable 
operating stress. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., July 17, 1970. 
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